Stability of standing matter waves in a trap 
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We discuss excited Bose-condensed states and find the criterion of dynamical stability of a kink-wise 
state, i.e., a standing matter wave with one nodal plane perpendicular to the axis of a cylindrical trap. 
The dynamical stability requires a strong radial confinement corresponding to the radial frequency 
larger than the mean-field interparticle interaction. We address the question of thermodynamic 
instability related to the presence of excitations with negative energy. 
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The discovery of Bose-Einstcin condensation (BEC) in 
trapped clouds of alkali atoms Q and extensive studies 
of Bose-condensed gases in the recent years have led to 
the observation of new macroscopic quantum phenom- 
ena, such as interference between two independently cre- 
ated condensates Q and reduction of the rates of inelastic 
processes (three-body recombination) in the presence of 
a condensate || . The success in these studies stimulates 
an interest in macroscopically excited Bose-condensed 
states, i.e., excited-state solutions of the Gross-Pitaevskii 
equation, where one can expect to observe novel signa- 
tures of BEC. A widely discussed example is the vortex 
state, well known in superfiuid liquid helium Q. 

Another option concerns excited Bose-condensed 
states which have a macroscopic wavefunction with nodal 
planes perpendicular to the symmetry axis of a trap. 
These states represent standing matter waves for which 
the trap serves as a cavity, and it will be worth studying 
in which aspects they are similar to light waves. An inter- 
esting idea concerns an atom laser for the generation of 
coherent matter waves which, due to the potential pres- 
ence of nodal planes in the condensate wavefunction, will 
be quite different from the matter waves out of a ground- 
state Bose condensate. The waves with one nodal plane 
(kinks or dark solitons), being the lowest energy phase 
slip, are of great interest in connection with the decay of 
persistent currents. Suggested ways of creating standing 
matter waves in a trap rely on adiabatic Raman transfer 
of particles from the ground to excited Bose-condensed 
state H or on selective population of trap levels by boson- 
ically enhanced spontaneous emission of optically excited 
atoms of an incoming beam Q . 

A principal question concerns the stability of excited 
Bose-condensed states with respect to the interparticle 
interaction. In this paper we consider standing matter 
waves with one nodal plane perpendicular to the axis of 
a cylindrical trap. For the axially Thomas-Fermi regime 
(axial frequency oj z is much smaller than the mean-field 
interaction) this state can be called "kink- wise" (see |7j 
and Fig.l), as the presence of the nodal plane makes a 
kink in the dependence of the condensate wavefunction 



^>o on the axial coordinate. A characteristic size of the 
kink is of order the correlation length, and the corre- 
sponding (axial) kinetic energy of the condensate is of 
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FIG. 1. Condensate wavefunction $0 of a kink- wise state 
in a cylindrical trap for the ratio of the mean-field inter- 
action at maximum density (nom) to the trap frequencies, 
nomU : hujp : huj z = 10 : 10 : 1, (U = 4irfi 2 a/m, a > is 
the scattering length, and m the atom mass). The axial (z) 
and radial (a;) coordinates are given in units of the correlation 

length I = %/ \J nomU, and (normalized to unity) ^0 in units 
of r 3/2 . In (b) the solid curve is ^0(0, z), and the dashed 
curve the first (axially) excited state of a harmonic oscillator. 
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order the mean-field interaction. Similarly to the case 
of vortices , the instability of the kink- wise state is 
related to the motion of the kink (core) with respect to 
the rest of the condensate. We analyze the spectrum of 
elementary excitations of this Bose-condensed state and 
find the criterion of dynamical stability, i.e., the stability 
of small-amplitude normal modes: In order to prevent 
the interaction-induced transfer of (axial) kinetic energy 
of the condensate to the radial degrees of freedom one 
should strongly confine the radial motion by making the 
radial frequency to p larger than the mean-field interparti- 
cle interaction. Under this condition the kink-wise state 
will be perfectly stable in the limit of zero temperature: 
The thermodynamic instability related to the presence of 
an excitation mode with negative energy will not lead to 
decay in the absence of dissipative processes. 

Our conclusion of how tightly one should confine the 
radial motion to achieve (quasi) ID dynamics of a conden- 
sate and observe dynamically stable kinks (dark solitons) 
is directly related to the problem of " engineering the di- 
mensionality of space" . Various experiments aiming for 
quasi-lD gases are currently being setup. 

To gain insight in the nature of the instability, we first 
consider a kink-wise Bose-condensed state in the absence 
of trapping field, i.e., the state with one nodal x, y plane 
in an otherwise spatially homogeneous condensate of den- 
sity no- We will use the chemical potential fj, = uqU and 
the correlation length / = h/^/mfi as units of energy and 
length, and no as unit of density. For a positive scattering 
length the Gross-Pitaevskii equation is reduced to 



1 d 2 * 

2 dz 2 



*o - *o = 



(1) 



and has a simple solution describing the kink in the de- 
pendence of \l/o on the z coordinate (see, e.g., Q|): 



= tanh z, 



(2) 



where z is counted from the nodal plane. The presence 
of the kink results in a large kinetic energy K (~ ji) per 
condensate particle at distances of order I from the nodal 
plane, which we below call kink-related kinetic energy: 



K(z) = 



1 T d 2 *o 



dz 2 



tanh 2 z 
cosh 2 z 



(3) 



Elementary excitations of the kink-wise condensate are 
characterized by the momentum k of transverse (x,y) 
motion and by the quantum number v of motion along 
the z axis. In terms of the operators bkv, V kv of the excita- 
tions the above-condensate part of the field operator of 
atoms can be represented as ^ fcl , exp(ikr)[ufcj, (£:)&£„ — 
Vkv{z)b_ ku ], and we obtain the following Bogolyubov- 
de Gennes equations for the excitation energies e^v and 
wavefunctions = utu ± Vkv'- 



e k uf T (z) = [h± + (A 2 /2)]/ ± (z), 



(4) 



where the operators h+ = —d 2 /2dz 2 
-d 2 /2dz 2 + 3^1 - 1 take the form 



1 and fi- 



ll- = - 



coth z d 
2 dz 

cosh 2 z d 
2 dz 



, o d 
tanh z— cothz 



dz 



1 



cosh z dz 



■ cosh z 



(5) 
(6) 



The kink-wise state has several excitation modes with 
zero energy. The ones for k — are the well-known 
fundamental modes of a ID kink state, following from 
decoupled equations h±f ± = 0. Those which do not 
exponentially grow at large z are = tanhz, /{" = 0; 
fi = 0' f-2 = cosh -2 z, and /g~ = ztanhz — 1, /£" = 0. 
The mode f% , which is even with respect to inversion of 
the z coordinate, is localized at distances of order I from 
the nodal plane. 

As an example, we demonstrate the instability of trans- 
verse normal modes in the absence of translational mo- 
tion of the nodal plane. We consider modes which in the 
limit k — > correspond to the localized zero-energy mode 
/ 2 . For transverse momenta k -C 1 the excitation wave- 
functions at distances z <C can be found as series of 
expansion in powers of A:. As the leading term in the ex- 
pansion for the function / _ is proportional to cosh -2 z, 
from Eqs.(^) one can find that the leading terms in the 
expansion for the function /+ should be /+ = const and 
/+ = /jj . Hence, this expansion takes the form 



/+ = C[l + Bf+(z) + k 2 g(z)}, 



(7) 



where C and B are constants. The equation for the func- 
tion g follows directly from Eqs. (jj)-©: 

2e\ = k 2 h-{2h + g(z) + Bf+(z) + 1) - k 2 cosh -2 z. (8) 

Using Eq.([j]) and performing the integration of Eq.(||), we 
obtain the relation {3e 2 k /k 2 + 1) cosh 2 z + 4 + 3/ 3 + (z) + 
6h + g(z) = 0. As g should not contain terms exponen- 
tially growing with z, we find the dispersion relation cor- 
responding to imaginary excitation energies: 



(9) 



The instability of transverse normal modes, following 
from Eq.(||), originates from the transfer of the (longi- 
tudinal) kink-related kinetic energy K of the condensate 
to these modes. As K ~ fj,, it can be transferred by the 
mean- field interaction to modes with small k. 

The demonstrated instability and Eq. (O) are similar 
to those in the case of "domain walls" [¥0[. In Fig. 2 
we present numerical results for Ime^. as a function of 
k. For small momenta it increases linearly with fc, in 
accordance with Eq. (|9|) , and reaches its maximum at k = 
l/v2- Further increase of k leads to decreasing Ime/. 
which becomes zero at k = 1. At this critical point we 
have one more zero-energy solution of Eqs.(^J): 



£fe = 0, / + oc 1/ coshz, / = 0. 



(10) 
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For k > 1 the energy of free transverse motion, fc 2 /2, ex- 
ceeds the kink-related kinetic energy K, and the normal 
modes are dynamically stable, with positive 




FIG. 2. Imaginary part of the excitation energy (in units 
of fj,) versus the transverse momentum k (in units of I' 1 ) for 
a kink- wise condensate in the absence of a trapping field. 

One can now understand the origin of dynamical insta- 
bility of a kink- wise Bose-condensed state in a cylindrical 
harmonic trap: The interparticle interaction can transfer 
the (axial) kink-related kinetic energy of the condensate 
to the radial degrees of freedom. In order to suppress 
this instability one has to significantly confine the radial 
motion. As the (axial) kinetic energy per particle in the 
axially Thomas-Fermi condensate is of order the mean- 
field interaction at maximum density, no m U, the radial 
frequency should be the same or larger. 

We have performed calculations for various ratia of the 
radial to axial frequency, u> p /u> z , and found the maximum 
value 7 C of the parameter 7 = no m U /%uj p , at which the 
kink-wise Bose-condensed state is still dynamically sta- 
ble, i.e., all excitation modes have real frequencies. If 
7 > 7c, there are excitations with imaginary frequencies, 
and the kink-wise condensate is dynamically unstable. 

We have solved the Gross-Pitaevskii equation 



2m 



A+^(LU 2 z z 2 +u;y) + U\^o\ 



*o = 0, 



(11) 



together with the Bogolyubov-de Gennes equations for 
the excitations, which we write in the form 



s2 r 



2m 



-A- 



/ ± + (lTl)C/|* | 2 / ± . (12) 



Eqs. ( |l2| ) give real e 2 which depend continuously on 7 
and the aspect ratio. In the range of 7 and w p /w z , where 
a given mode v is dynamically unstable, e 2 v < and the 
energy e v is purely imaginary. In the region of dynamical 
stability e v is purely real (el > 0), and, hence, at the 
border between the two regions we have s u = 0. 

At the critical point 7 = 7 C all excitation energies e v 
are real, and one of the excitations has zero energy. This 



is just the mode which for 7 > j c becomes dynamically 
unstable. Similarly to the mode of Eq.(|l0|) in the ab- 
sence of trapping field, this mode is even with respect to 
inversion of the z coordinate. The function / _ = 0, and 
/ + follows directly from Eqs.(|l2]): 



(-A + A* /*o)/ + = 0. 



(13) 



Eq.(|13[) is the Schrodinger equation for the motion of a 
particle (with zero energy) in a cylindrically symmetric 
potential V — h 2 A 1 i'o/2m 1 ^Q. The potential V depends 
on 7 and the aspect ratio. Thus, for a given ratio to p /u) z 
one finds the critical value j c by selecting the param- 
eter 7 such that there is an even (non-zero) solution of 
Eq. ( |l3| ) , remaining finite at the origin and tending to zero 
at infinity. This was checked numerically on the basis of 
Eqs. (pr|)-(p^|) for a wide range of 7 and the aspect ratio. 

As follows from our calculations, j c is minimal for exci- 
tations with the projection of the orbital angular momen- 
tum on the symmetry axis, M = 1. The dependence of 
7 C on the aspect ratio is presented in Fig. 3. For uj p < lu z 
even an arbitrary small interparticle interaction leads to 
instability, as the axial "kink-related" energy per parti- 
cle in the condensate (Tiuj z ) can be always transferred to 
the radial mode with M = 1 , which by itself has energy 
Tiujp. For u) p > lu z the critical value 7 C increases with the 
ratio uj p /oj z and reaches 7 C ~ 2.4 for uj p 3> w z . We also 
found that the decay of dynamically unstable kink states 
is accompanied by the undulation of the nodal plane and 
the formation of vortex-antivortex pairs, similar to the 
decay of dark optical solitons (ll) . 




FIG. 3. Critical parameter 7 C versus the aspect ratio for 
a kink- wise condensate in a cylindrical trap. 

The criterion of dynamical stability of a kink- wise con- 
densate, 7 < 7c, can be satisfied in the conditions of cur- 
rent BEC experiments. For a rubidium condensate in a 
cylindrical trap with lo p ~ 200 Hz 3> lo z it requires the 
maximum density no m < 10 cm~ 3 . 

Although for 7 < 7 C the kink-wise condensate is dy- 
namically stable, there is a thermodynamic instability 
related to the presence of an excitation with negative en- 
ergy. For a very strong radial confinement of the axially 
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Thomas- Fermi kink- wise condensate (Tbui p 3> na m U ^> 
huj z ; 7 < 7 C ) we calculate a negative excitation en- 
ergy close to e* = —Tiu z /\/2 characteristic for the 1-D 
Thomas-Fermi kink-wise condensate in a harmonic trap. 

In the 1-D case we calculate the negative excita- 
tion energy analytically, by solving the Bogolyubov-de 
Gennes equations at distances z from the origin, much 
smaller than the Thomas-Fermi size of the condensate 
R = {2^/muj 2 ) 1 / 2 . We represent and the excitation 
wavefunctions as series of expansion in powers of small 
parameter £ = tku z /fi. Then, in the same dimension- 
less units as in the absence of trapping field, the Gross- 
Pitaevskii equation is given by Eq.(|IJ), with an extra term 
C, 2 z 2 ^Iq/2 in the left-hand side. Confining ourselves to 
the expansion up to £ 2 , we obtain 



*o = tanhz + C 2 ??(z), 



(14) 



where the function r/(z) is determined by the equation 

h-T](z) + 2 /2)tanhz = 0, (15) 

and is not given because of its complexity. For \z\ ^> 1 
we have r\ = — signal + 2z 2 )/8, and Eq. (|l4|) recovers the 
Thomas-Fermi result at \z\ -C R/l. The Bogolyubov-de 
Gennes equations take the form 

s»f T = h ± f ± + C 2 [z 2 /2 + (4 T 2Hz) tanh z]f± (16) 

Just as in the absence of trapping field, we consider a 
mode for which the leading term in the expansion for 
the function / _ is proportional to 1/ cosh 2 z. For finding 
the excitation energy it is sufficient to keep the terms 
independent of £ and proportional to C, 2 in the expansion 
for the function / + . Then, similarly to Eq.(^), we obtain 
/+ oc [1 + ( 2 G(z)]. The equation for G(z) follows from 
Eqs. (Il6]), and by using Eqs. (p|),(p"5|) is transformed to 



1 2 

, 9 d 
2 cosh z — 
dz 



h-[h+G(z) + z 2 
tanh z d 



cosh z dz 



(cosh zrj(z)) 



(17) 



d 2 G/dz 



Integration of Eq.(|l7|) gives at large \z\ the relation 



(e*/C - l/2)cosh J z - 1 and, as G should 



not contain exponentially growing terms, we obtain e 2 = 
C 2 /2. The normalization condition J dzf + f~ = 1 allows 
then to conclude that the excitation energy is negative 
and, hence, equal to —%u> z /\/2. The frequency u z /y2 
for the oscillations of the kink in a 1-D Thomas-Fermi 
condensate has also been found in the recent work |fjj] 
from the equation of motion for the kink. 

The excitation spectrum for e v > follows from the 
solution of the Bogolubov-de Gennes equations at large 
\z\, where the Thomas-Fermi shape of |^o| is not influ- 
enced by the kink. Then, along the lines of the theory 
for the 3D case [na Of , we obtain a discrete spectrum 



Finally, we analyze the influence of the excitation with 
negative energy on the stability of the kink- wise conden- 
sate. Beyond the Bogolubov-de Gennes approach, there 
is a small coupling of this excitation to the excitations 
with positive energies. But at temperatures T —> there 
will be no real decay processes. Those require a simulta- 
neous creation of excitations with positive and negative 
energies, with the total excitation energy equal to zero. 
Due to the structure of the discrete spectrum for e v > 0, 
this conservation of energy can not be satisfied while cre- 
ating a moderate number of excitations with the above 
found negative energy —hio z j \[2. 

Thus, under the condition of dynamical stability the 
kink-wise Bose-condensed state is perfectly stable at 
T — > 0. The decay mechanism in the presence of a ther- 
mal cloud is in some sense similar to that of temperature- 
dependent damping of excitations in trapped Bose- 
condensed gases (see |l6|]) and originates from the scat- 
tering of thermal particles on the kink. Accordingly, the 
decay time can be made large by decreasing temperature 
well below the value of the mean-field interparticle inter- 
action. A detailed analysis of dissipative dynamics of a 
kink state at finite T requires a separate investigation. 
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e v = huo z yJ v{v + l)/2, where v is a positive integer. 
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